Genome-wide Analysis Identifies Nuclear Factor 1C as a Novel Transcription Factor and Potential Therapeutic Target in Small Cell Lung Cancer.

BACKGROUND: Recent insights regarding mechanisms mediating stemness, heterogeneity, and metastatic potential of lung cancers have yet to be fully translated to effective regimens for the treatment of these malignancies. This study sought to identify novel targets for lung cancer therapy. METHODS: Transcriptomes and DNA methylomes of 14 SCLC and 10 NSCLC lines were compared to normal human small airway epithelial cells (SAEC) and induced pluripotent stem cell (iPSC) clones derived from SAEC. SCLC lines, lung iPSC (Lu-iPSC), and SAEC were further evaluated by DNase I hypersensitivity (DHS-seq). Changes in chromatin accessibility and depths of transcription factor (TF) footprints were quantified using Bivariate analysis of Genomic Footprint. Standard techniques were used to examine growth and tumorigenencity as well as changes in transcriptomes and glucose metabolism of SCLC cells following Nuclear Factor 1C (NFIC) knockdown, and to examine NFIC expression in SCLC cells following exposure to BET inhibitors. RESULTS: Significant commonality of transcriptomes and DNA methylomes was observed between Lu-iPSC and SCLC; however, this analysis was uninformative regarding pathways unique to lung cancer. Linking results of DNase-seq to RNA-seq enabled identification of networks not previously associated with SCLC. When combined with footprint depth, NFIC, a transcription factor not previously associated with SCLC, had the highest score of occupancy at open chromatin sites. Knockdown of NFIC impaired glucose metabolism, decreased stemness, and inhibited growth of SCLC cells in-vitro and in-vivo. ChIP-seq analysis identified numerous sites occupied by Bromodomain-containing protein 4 (BRD4) in the NFIC promoter region. Knock-down of BRD4 or treatment with Bromodomain and extra-terminal domain (BET) inhibitors (BETi) markedly reduced NFIC expression in SCLC cells and SCLC PDX models. Approximately 8% of genes downregulated by BETi treatment were repressed by NFIC knockdown in SCLC, while 34% of genes repressed following NFIC knockdown were also downregulated in SCLC cells following BETi treatment. CONCLUSIONS: NFIC is a key TF and possible mediator of transcriptional regulation by BET family proteins in SCLC. Our findings highlight the potential of genome-wide chromatin accessibility analysis for elucidating mechanisms of pulmonary carcinogenesis and identifying novel targets for lung cancer therapy.

Novel biosensor for high-throughput detection of progesterone receptor-interacting endocrine disruptors.

Progesterone receptor (PR)-interacting compounds in the environment are associated with serious health hazards. However, methods for their detection in environmental samples are cumbersome. We report a sensitive activity-based biosensor for rapid and reliable screening of progesterone receptor (PR)-interacting endocrine disrupting chemicals (EDCs). The biosensor is a cell line which expresses nuclear mCherry-NF1 and a green fluorescent protein (GFP)-tagged chimera of glucocorticoid receptor (GR) N terminus fused to the ligand binding domain (LBD) of PR (GFP-GR-PR). As this LBD is shared by the PRA and PRB, the biosensor reports on the activation of both PR isoforms. This GFP-GR-PR chimera is cytoplasmic in the absence of hormone and translocates rapidly to the nucleus in response to PR agonists or antagonists in concentration- and time-dependent manner. In live cells, presence of nuclear NF1 label eliminates cell fixation and nuclear staining resulting in efficient screening. The assay can be used in screens for novel PR ligands and PR-interacting contaminants in environmental samples. A limited screen of river water samples indicated a widespread, low-level contamination with PR-interacting contaminants in all tested samples.

The mineralocorticoid receptor forms higher order oligomers upon DNA binding.

The prevailing model of steroid hormone nuclear receptor function assumes ligand-induced homodimer formation followed by binding to DNA hormone response elements (HREs). This model has been challenged by evidence showing that the glucocorticoid receptor (GR) forms tetramers upon ligand and DNA binding, which then drive receptor-mediated gene transactivation and transrepression. GR and the closely-related mineralocorticoid receptors (MR) interact to transduce corticosteroid hormone signaling, but whether they share the same quaternary arrangement is unknown. Here, we used a fluorescence imaging technique, Number & Brightness, to study oligomerization in a cell system allowing real-time analysis of receptor-DNA interactions. Agonist-bound MR forms tetramers in the nucleoplasm and higher order oligomers upon binding to HREs. Antagonists form intermediate-size quaternary arrangements, suggesting that large oligomers are essential for function. Divergence between MR and GR quaternary structure is driven by different functionality of known and new multimerization interfaces, which does not preclude formation of heteromers. Thus, influencing oligomerization may be important to selectively modulate corticosteroid signaling.

Transcription Factor Dynamics: One Molecule at a Time.

Cells must tightly regulate their gene expression programs and yet rapidly respond to acute biochemical and biophysical cues within their environment. This information is transmitted to the nucleus through various signaling cascades, culminating in the activation or repression of target genes. Transcription factors (TFs) are key mediators of these signals, binding to specific regulatory elements within chromatin. While live-cell imaging has conclusively proven that TF-chromatin interactions are highly dynamic, how such transient interactions can have long-term impacts on developmental trajectories and disease progression is still largely unclear. In this review, we summarize our current understanding of the dynamic nature of TF functions, starting with a historical overview of early live-cell experiments. We highlight key factors that govern TF dynamics and how TF dynamics, in turn, affect downstream transcriptional bursting. Finally, we conclude with open challenges and emerging technologies that will further our understanding of transcriptional regulation.

Long-range gene regulation in hormone-dependent cancer.

The human genome is organized into multiple structural layers, ranging from chromosome territories to progressively smaller substructures, such as topologically associating domains (TADs) and chromatin loops. These substructures, collectively referred to as long-range chromatin interactions (LRIs), have a significant role in regulating gene expression. TADs are regions of the genome that harbour groups of genes and regulatory elements that frequently interact with each other and are insulated from other regions, thereby preventing widespread uncontrolled DNA contacts. Chromatin loops formed within TADs through enhancer and promoter interactions are elastic, allowing transcriptional heterogeneity and stochasticity. Over the past decade, it has become evident that the 3D genome structure, also referred to as the chromatin architecture, is central to many transcriptional cellular decisions. In this Review, we delve into the intricate relationship between steroid receptors and LRIs, discussing how steroid receptors interact with and modulate these chromatin interactions. Genetic alterations in the many processes involved in organizing the nuclear architecture are often associated with the development of hormone-dependent cancers. A better understanding of the interplay between architectural proteins and hormone regulatory networks can ultimately be exploited to develop improved approaches for cancer treatment.

SMYD3 represses tumor-intrinsic interferon response in HPV-negative squamous cell carcinoma of the head and neck.

Cancers often display immune escape, but the mechanisms are incompletely understood. Herein, we identify SMYD3 as a mediator of immune escape in human papilloma virus (HPV)-negative head and neck squamous cell carcinoma (HNSCC), an aggressive disease with poor response to immunotherapy with pembrolizumab. SMYD3 depletion induces upregulation of multiple type I interferon (IFN) response and antigen presentation machinery genes in HNSCC cells. Mechanistically, SMYD3 binds to and regulates the transcription of UHRF1, encoding for a reader of H3K9me3, which binds to H3K9me3-enriched promoters of key immune-related genes, recruits DNMT1, and silences their expression. SMYD3 further maintains the repression of immune-related genes through intragenic deposition of H4K20me3. In vivo, Smyd3 depletion induces influx of CD8+ T cells and increases sensitivity to anti-programmed death 1 (PD-1) therapy. SMYD3 overexpression is associated with decreased CD8 T cell infiltration and poor response to neoadjuvant pembrolizumab. These data support combining SMYD3 depletion strategies with checkpoint blockade to overcome anti-PD-1 resistance in HPV-negative HNSCC.

Dynamic switching of transcriptional regulators between two distinct low-mobility chromatin states.

How chromatin dynamics relate to transcriptional activity remains poorly understood. Using single-molecule tracking, coupled with machine learning, we show that histone H2B and multiple chromatin-bound transcriptional regulators display two distinct low-mobility states. Ligand activation results in a marked increase in the propensity of steroid receptors to bind in the lowest-mobility state. Mutational analysis revealed that interactions with chromatin in the lowest-mobility state require an intact DNA binding domain and oligomerization domains. These states are not spatially separated as previously believed, but individual H2B and bound-TF molecules can dynamically switch between them on time scales of seconds. Single bound-TF molecules with different mobilities exhibit different dwell time distributions, suggesting that the mobility of TFs is intimately coupled with their binding dynamics. Together, our results identify two unique and distinct low-mobility states that appear to represent common pathways for transcription activation in mammalian cells.

The mineralocorticoid receptor forms higher order oligomers upon DNA binding.

The prevailing model of steroid hormone nuclear receptor function assumes ligand-induced homodimer formation followed by binding to DNA hormone response elements (HREs). This model has been challenged by evidence showing that the glucocorticoid receptor (GR) forms tetramers upon ligand and DNA binding, which then drive receptor-mediated gene transactivation and transrepression. GR and the closely-related mineralocorticoid receptors (MR) interact to transduce corticosteroid hormone signaling, but whether they share the same quaternary arrangement is unknown. Here, we used a fluorescence imaging technique, Number & Brightness, to study oligomerization in a cell system allowing real-time analysis of receptor-DNA interactions. Agonist-bound MR forms tetramers in the nucleoplasm and higher order oligomers upon binding to HREs. Antagonists form intermediate quaternary arrangements, suggesting that large oligomers are essential for function. Divergence between MR and GR quaternary structure is driven by different functionality of known and new multimerization interfaces, which does not preclude formation of heteromers. Thus, influencing oligomerization may be important to selectively modulate corticosteroid signaling.

The Glucocorticoid Receptor is Required for Efficient Aldosterone-Induced Transcription by the Mineralocorticoid Receptor.

The glucocorticoid and mineralocorticoid receptors (GR and MR, respectively) have distinct, yet overlapping physiological and pathophysiological functions. There are indications that both receptors interact functionally and physically, but the precise role of this interdependence is poorly understood. Here, we analyzed the impact of GR co-expression on MR genome-wide chromatin binding and transcriptional responses to aldosterone and glucocorticoids, both physiological ligands of this receptor. Our data show that GR co-expression alters MR genome-wide binding to consensus DNA sequences in a locus- and ligand-specific way. MR binding to consensus DNA sequences is affected by GR. Transcriptional responses of MR in the absence of GR are weak and show poor correlation with chromatin binding. In contrast, co-expression of GR potentiates MR-mediated transcription, particularly in response to aldosterone. Finally, single-molecule tracking of MR suggests that the presence of GR contributes to productive binding of MR/aldosterone complexes to chromatin. Together, our data indicate that co-expression of GR potentiates aldosterone-mediated MR transcriptional activity, even in the absence of glucocorticoids.

HIV vaccine candidate efficacy in female macaques mediated by cAMP-dependent efferocytosis and V2-specific ADCC.

The development of an effective vaccine to protect against HIV acquisition will be greatly bolstered by in-depth understanding of the innate and adaptive responses to vaccination. We report here that the efficacy of DNA/ALVAC/gp120/alum vaccines, based on V2-specific antibodies mediating apoptosis of infected cells (V2-ADCC), is complemented by efferocytosis, a cyclic AMP (cAMP)-dependent antiphlogistic engulfment of apoptotic cells by CD14+ monocytes. Central to vaccine efficacy is the engagement of the CCL2/CCR2 axis and tolerogenic dendritic cells producing IL-10 (DC-10). Epigenetic reprogramming in CD14+ cells of the cyclic AMP/CREB pathway and increased systemic levels of miRNA-139-5p, a negative regulator of expression of the cAMP-specific phosphodiesterase PDE4D, correlated with vaccine efficacy. These data posit that efferocytosis, through the prompt and effective removal of apoptotic infected cells, contributes to vaccine efficacy by decreasing inflammation and maintaining tissue homeostasis.

Juxtaposition of intensive agriculture, vulnerable aquifers, and mixed chemical/microbial exposures in private-well tapwater in northeast Iowa.

In the United States and globally, contaminant exposure in unregulated private-well point-of-use tapwater (TW) is a recognized public-health data gap and an obstacle to both risk-management and homeowner decision making. To help address the lack of data on broad contaminant exposures in private-well TW from hydrologically-vulnerable (alluvial, karst) aquifers in agriculturally-intensive landscapes, samples were collected in 2018-2019 from 47 northeast Iowa farms and analyzed for 35 inorganics, 437 unique organics, 5 in vitro bioassays, and 11 microbial assays. Twenty-six inorganics and 51 organics, dominated by pesticides and related transformation products (35 herbicide-, 5 insecticide-, and 2 fungicide-related), were observed in TW. Heterotrophic bacteria detections were near ubiquitous (94 % of the samples), with detection of total coliform bacteria in 28 % of the samples and growth on at least one putative-pathogen selective media across all TW samples. Health-based hazard index screening levels were exceeded frequently in private-well TW and attributed primarily to inorganics (nitrate, uranium). Results support incorporation of residential treatment systems to protect against contaminant exposure and the need for increased monitoring of rural private-well homes. Continued assessment of unmonitored and unregulated private-supply TW is needed to model contaminant exposures and human-health risks.

Extrachromosomal DNA Amplification Contributes to Small Cell Lung Cancer Heterogeneity and Is Associated with Worse Outcomes.

Small-cell lung cancer (SCLC) is an aggressive neuroendocrine lung cancer. Oncogenic MYC amplifications drive SCLC heterogeneity, but the genetic mechanisms of MYC amplification and phenotypic plasticity, characterized by neuroendocrine and nonneuroendocrine cell states, are not known. Here, we integrate whole-genome sequencing, long-range optical mapping, single-cell DNA sequencing, and fluorescence in situ hybridization to find extrachromosomal DNA (ecDNA) as a primary source of SCLC oncogene amplifications and driver fusions. ecDNAs bring to proximity enhancer elements and oncogenes, creating SCLC transcription-amplifying units, driving exceptionally high MYC gene dosage. We demonstrate that cell-free nucleosome profiling can noninvasively detect ecDNA amplifications in plasma, facilitating its genome-wide interrogation in SCLC and other cancers. Altogether, our work provides the first comprehensive map of SCLC ecDNA and describes a new mechanism that governs MYC-driven SCLC heterogeneity. ecDNA-enabled transcriptional flexibility may explain the significantly worse survival outcomes of SCLC harboring complex ecDNA amplifications. SIGNIFICANCE: MYC drives SCLC progression, but the genetic basis of MYC-driven SCLC evolution is unknown. Using SCLC as a paradigm, we report how ecDNA amplifications function as MYC-amplifying units, fostering tumor plasticity and a high degree of tumor heterogeneity. This article is highlighted in the In This Issue feature, p. 799.

The multivalency of the glucocorticoid receptor ligand-binding domain explains its manifold physiological activities.

The glucocorticoid receptor (GR) is a ubiquitously expressed transcription factor that controls metabolic and homeostatic processes essential for life. Although numerous crystal structures of the GR ligand-binding domain (GR-LBD) have been reported, the functional oligomeric state of the full-length receptor, which is essential for its transcriptional activity, remains disputed. Here we present five new crystal structures of agonist-bound GR-LBD, along with a thorough analysis of previous structural work. We identify four distinct homodimerization interfaces on the GR-LBD surface, which can associate into 20 topologically different homodimers. Biologically relevant homodimers were identified by studying a battery of GR point mutants including crosslinking assays in solution, quantitative fluorescence microscopy in living cells, and transcriptomic analyses. Our results highlight the relevance of non-canonical dimerization modes for GR, especially of contacts made by loop L1-3 residues such as Tyr545. Our work illustrates the unique flexibility of GR's LBD and suggests different dimeric conformations within cells. In addition, we unveil pathophysiologically relevant quaternary assemblies of the receptor with important implications for glucocorticoid action and drug design.

Identification of a novel GR-ARID1a-P53BP1 protein complex involved in DNA damage repair and cell cycle regulation.

ARID1a (BAF250), a component of human SWI/SNF chromatin remodeling complexes, is frequently mutated across numerous cancers, and its loss of function has been putatively linked to glucocorticoid resistance. Here, we interrogate the impact of siRNA knockdown of ARID1a compared to a functional interference approach in the HeLa human cervical cancer cell line. We report that ARID1a knockdown resulted in a significant global decrease in chromatin accessibility in ATAC-Seq analysis, as well as affecting a subset of genome-wide GR binding sites determined by analyzing GR ChIP-Seq data. Interestingly, the specific effects on gene expression were limited to a relatively small subset of glucocorticoid-regulated genes, notably those involved in cell cycle regulation and DNA repair. The vast majority of glucocorticoid-regulated genes were largely unaffected by ARID1a knockdown or functional interference, consistent with a more specific role for ARID1a in glucocorticoid function than previously speculated. Using liquid chromatography-mass spectrometry, we have identified a chromatin-associated protein complex comprising GR, ARID1a, and several DNA damage repair proteins including P53 binding protein 1 (P53BP1), Poly(ADP-Ribose) Polymerase 1 (PARP1), DNA damage-binding protein 1 (DDB1), DNA mismatch repair protein MSH6 and splicing factor proline and glutamine-rich protein (SFPQ), as well as the histone acetyltransferase KAT7, an epigenetic regulator of steroid-dependent transcription, DNA damage repair and cell cycle regulation. Not only was this protein complex ablated with both ARID1a knockdown and functional interference, but spontaneously arising DNA damage was also found to accumulate in a manner consistent with impaired DNA damage repair mechanisms. Recovery from dexamethasone-dependent cell cycle arrest was also significantly impaired. Taken together, our data demonstrate that although glucocorticoids can still promote cell cycle arrest in the absence of ARID1a, the purpose of this arrest to allow time for DNA damage repair is hindered.

Transcription factors TCF-1 and GATA3 are key factors for the epigenetic priming of early innate lymphoid progenitors toward distinct cell fates.

The differentiation of innate lymphoid cells (ILCs) from hematopoietic stem cells needs to go through several multipotent progenitor stages. However, it remains unclear whether the fates of multipotent progenitors are predefined by epigenetic states. Here, we report the identification of distinct accessible chromatin regions in all lymphoid progenitors (ALPs), EILPs, and ILC precursors (ILCPs). Single-cell MNase-seq analyses revealed that EILPs contained distinct subpopulations epigenetically primed toward either dendritic cell lineages or ILC lineages. We found that TCF-1 and GATA3 co-bound to the lineage-defining sites for ILCs (LDS-Is), whereas PU.1 binding was enriched in the LDSs for alternative dendritic cells (LDS-As). TCF-1 and GATA3 were indispensable for the epigenetic priming of LDSs at the EILP stage. Our results suggest that the multipotency of progenitor cells is defined by the existence of a heterogeneous population of cells epigenetically primed for distinct downstream lineages, which are regulated by key transcription factors.

The glucocorticoid receptor associates with the cohesin loader NIPBL to promote long-range gene regulation.

The cohesin complex is central to chromatin looping, but mechanisms by which these long-range chromatin interactions are formed and persist remain unclear. We demonstrate that interactions between a transcription factor (TF) and the cohesin loader NIPBL regulate enhancer-dependent gene activity. Using mass spectrometry, genome mapping, and single-molecule tracking methods, we demonstrate that the glucocorticoid (GC) receptor (GR) interacts with NIPBL and the cohesin complex at the chromatin level, promoting loop extrusion and long-range gene regulation. Real-time single-molecule experiments show that loss of cohesin markedly diminishes the concentration of TF molecules at specific nuclear confinement sites, increasing TF local concentration and promoting gene regulation. Last, patient-derived acute myeloid leukemia cells harboring cohesin mutations exhibit a reduced response to GCs, suggesting that the GR-NIPBL-cohesin interaction is defective in these patients, resulting in poor response to GC treatment.

Nuclear pore protein NUP210 depletion suppresses metastasis through heterochromatin-mediated disruption of tumor cell mechanical response.

Mechanical signals from the extracellular microenvironment have been implicated in tumor and metastatic progression. Here, we identify nucleoporin NUP210 as a metastasis susceptibility gene for human estrogen receptor positive (ER+) breast cancer and a cellular mechanosensor. Nup210 depletion suppresses lung metastasis in mouse models of breast cancer. Mechanistically, NUP210 interacts with LINC complex protein SUN2 which connects the nucleus to the cytoskeleton. In addition, the NUP210/SUN2 complex interacts with chromatin via the short isoform of BRD4 and histone H3.1/H3.2 at the nuclear periphery. In Nup210 knockout cells, mechanosensitive genes accumulate H3K27me3 heterochromatin modification, mediated by the polycomb repressive complex 2 and differentially reposition within the nucleus. Transcriptional repression in Nup210 knockout cells results in defective mechanotransduction and focal adhesion necessary for their metastatic capacity. Our study provides an important role of nuclear pore protein in cellular mechanosensation and metastasis.

Activation of Crtc2/Creb1 in skeletal muscle enhances weight loss during intermittent fasting.

The Creb-Regulated Transcriptional Coactivator (Crtc) family of transcriptional coregulators drive Creb1-mediated transcription effects on metabolism in many tissues, but the in vivo effects of Crtc2/Creb1 transcription on skeletal muscle metabolism are not known. Skeletal muscle-specific overexpression of Crtc2 (Crtc2 mice) induced greater mitochondrial activity, metabolic flux capacity for both carbohydrates and fats, improved glucose tolerance and insulin sensitivity, and increased oxidative capacity, supported by upregulation of key metabolic genes. Crtc2 overexpression led to greater weight loss during alternate day fasting (ADF), selective loss of fat rather than lean mass, maintenance of higher energy expenditure during the fast and reduced binge-eating during the feeding period. ADF downregulated most of the mitochondrial electron transport genes, and other regulators of mitochondrial function, that were substantially reversed by Crtc2-driven transcription. Glucocorticoids acted with AMPK to drive atrophy and mitophagy, which was reversed by Crtc2/Creb1 signaling. Crtc2/Creb1-mediated signaling coordinates metabolic adaptations in skeletal muscle that explain how Crtc2/Creb1 regulates metabolism and weight loss.

Corticosterone pattern-dependent glucocorticoid receptor binding and transcriptional regulation within the liver.

Ultradian glucocorticoid rhythms are highly conserved across mammalian species, however, their functional significance is not yet fully understood. Here we demonstrate that pulsatile corticosterone replacement in adrenalectomised rats induces a dynamic pattern of glucocorticoid receptor (GR) binding at ~3,000 genomic sites in liver at the pulse peak, subsequently not found during the pulse nadir. In contrast, constant corticosterone replacement induced prolonged binding at the majority of these sites. Additionally, each pattern further induced markedly different transcriptional responses. During pulsatile treatment, intragenic occupancy by active RNA polymerase II exhibited pulsatile dynamics with transient changes in enrichment, either decreased or increased depending on the gene, which mostly returned to baseline during the inter-pulse interval. In contrast, constant corticosterone exposure induced prolonged effects on RNA polymerase II occupancy at the majority of gene targets, thus acting as a sustained regulatory signal for both transactivation and repression of glucocorticoid target genes. The nett effect of these differences were consequently seen in the liver transcriptome as RNA-seq analysis indicated that despite the same overall amount of corticosterone infused, twice the number of transcripts were regulated by constant corticosterone infusion, when compared to pulsatile. Target genes that were found to be differentially regulated in a pattern-dependent manner were enriched in functional pathways including carbohydrate, cholesterol, glucose and fat metabolism as well as inflammation, suggesting a functional role for dysregulated glucocorticoid rhythms in the development of metabolic dysfunction.